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ABSTRACT 

We present transit observations of the WASP-2 exoplanet system by the Apache Point Survey 
of Transit Lightcurves of Exoplanet s (APOSTLE) program. Model fitting to these data allows us 
to improve measurements of the hot- Jupiter exoplanet WASP-2b and its orbital parameters by a 
factor of ~ 2 over prior studies; we do not find evidence for transit depth variations. We do find 
reduced % 2 values greater than 1.0 in the observed minus computed transit times. A sinusoidal fit 
to the residuals yields a timing semi-amplitude of 32 seconds and a period of 389 days. However, 
random rearrangements of the data provide similar quality fits, and we cannot with certainty ascribe 
the timing variations to mutual exoplanet interactions. This inconclusive result is consistent with 
the lack of incontrovertible transit timing variations (TTVs) observed in other hot-Jupiter systems. 
This outcome emphasizes that unique recognition of TTVs requires dense sampling of the libration 
cycle (e.g. continuous observations from space-based platforms). However, even in systems observed 
with the Kepler spacecraft, there is a noted lack of transiting companions and TTVs in hot— Jupiter 
systems. This result is more meaningful, and indicates that hot- Jupiter systems, while they are easily 
observable from the ground, do not appear to be currently configured in a manner favorable to the 
detection of TTVs. The future of ground-based TTV studies may reside in resolving secular trends, 
and/or implementation at extreme quality observing sites to minimize atmospheric red noise. 

Subject headings: eclipses, stars: planetary systems, planets and satellites: fundamental parameters, 
individual: WASP-2b 



1. INTRODUCTION 

The transit technique is a highly efficient means of 
searching for exoplanetary systems. The detection of 
transiting systems requires a fortuitous observing geom- 
etry during the experiment - the exoplanet must be ob- 
served to traverse its host's stellar disk, occurring for 
~ 10% of viewing angles for hot Jupiters, but only 
0.4% for Earth analogues - and a control of experi- 
mental systematics at or below the level of the tran- 
sit depth. For hot- Jupiter systems around Solar-type 
stars, this may be as high as 1% (10 4 parts per mil- 
lion, ppm) of the out-of-transit depth; an Earth ana- 
logue in the same system would cause a transit depth 
of only 85 ppm. Ground-based observations have been 
able t o achieve per-exposure precisions o f down to 211 
ppm (iTregloan-Reed fc Southworthl 120121) to 250 ppm 
( Gillil and et al.l 119931 ) for the brightest objects. More 
commonly, in transit follow-up efforts where measure- 
ments of faint stars in the field is not a priority, rela- 
tive photometry at the 300-500 ppm level is achieved 
through defocusing and precision tracking, which mini- 
mize the sampling of the flat-field function, and allow 
the observer to accumulate more photons before sat- 
uration (e.g. iSouthworth et al l 120091 iWinn et al.l 120091 : 
iGillon et al1l20l2t iLendl et alJ^Oll T Methods to model 
the non-random ( "red" ) noise in the data are also shown 
to impr ove the accuracy of ph otometric transit measure- 
ments (|Carter fc Winn] I2009H . However, relative preci- 
sion in ground-based data is ultimately limited by atmo- 
spheric decoherence between the target and comparison 
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stars and across observing epochs, with le ading terms 
inclu ding the structure function of clouds (jlvezic et al.l 
I2007T ). the time-rate of change of aerosols a nd water va- 
por in the atmosphe re dStubbs et all 120071) . and atmo- 
spheric scintillation (|Young et al.l I1991D . Space-based 
obser vations by the Kepler spacecraft (jBorucki et al.l 
12010( 1 set the gold-standard for relative photometry, 
reaching root-mean-square syste matic variations of 2 
ppm on timescales of several hours (|Gilhland et al J 20 111 ) . 
As recognized by lAgol et all (|2005l ) and 



iHolman fc Murravl (|2005l ). the times of transits in 
multi-planet systems may not be exactly periodic due 
to mutual gravitational interactions of the planets. 
Multiple studies have been undertaken to search for 
these transit timing variations (TTVs) in known tran- 
siting exoplanet systems, using high-precision follow-up 
observations. The vast majority of follow-up has 
been undertaken on systems originally discovered from 
the ground, which have been heavily biased towards 
hot- Jupiter systems. The follow-up sampling of these 
transits is often irregular, due to weather, daytime, 
and seasonal effects. This makes detection of the 
libration of transit times, expected to occur over the 
timescales of months to years, dif ficult to recognize 
because the signal is undersampled. iBoue et all (|2012l ) 
outline the difficulties in resolving TTVs in Jovian 
systems: they are mostly relevant for systems near 
(but not exactly at) mean-motion resonance, and even 
detected signals may yield degenerate solutions for 
the mass of the perturber. No unambiguous TTVs 
have yet bee n discovered in gro und-based follow-up 
(see however iTinglev et al.1 120111 ). In contrast, the 
Kepler spacecraft follows-up its own discoveries through 
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continuous lightcurve coverage. Several multi-planet 
systems have bee n discovered th r ough Kepler TTVs (e.g. 
Ford et al.l |20H iSteffen et all l2012at iFabrvckv etall 



APOSTLE Observations of WASP2 



2012ft . which have proven to be a powerful verification 



and mass measureme nt technique for planet candidates 
(|Cochran et al.ll2011h . 

2. APOSTLE PROGRAM 

The Apache Point S urvey of Transit Lightcur ves of Ex- 
oplanets (APOSTLE: iKundurthv et al.ll2013bh program 
was initiated as a systematic study of known transiting 
cxoplan et systems on the ARC 3.5m telescope + Agile 
imager (M ukadam et al.l 1201 ID . The large aperture of 
the system and frame-transfer capabilities of Agile al- 
low us to obtain high-precision (500 ppm RMS) relative 
photometry between R = 10.8 th and R = 10 . 8*^ mag- 
nitude stars (XO-2; IKundurthv et al.l [2013a'). decreas- 
ing to 700 ppm at R = 12.2 vs. R = 12.9 (TrES-3; 
IKundurthv etall l2013bl ). and 1000 ppm at R = 13.8 
vs. R = 13.6 (GJ 1214; IKundurthv et alJl201lh . Im- 
portantly, these observations happen at 100% duty cycle 
due to Agile's frame-transfer capabilities. We make use 
of the Mult iTrans it Quick modeling program described 
in IKundurthv et all (|2013bf ) . which uses a Markov Chain 
Monte Carlo (MCMC) analyzer alongside a lightcurve 
parameterization that minimizes degeneracies between 
fitted parameters, ensuring that the MCMC proceeds ef- 
ficiently and faithfully samples parameter space. 

This pa per describes APOST LE observations of the 
WASP-2 (|Cameron et all 120071 ) system. The host star 
WASP-2A is a R = 11.3, spectral type Kl dwarf, 
with an effective temperature of T e ff = 5110 ± 60 in- 
ferred from optical and infrared colors (jMaxted et al.1 
l201lh . T eff = 5150 ± 80K using photos pheric fitting 
of spectroscopic data (jTriaud et al.l 12010ft. and m etal- 
licity of [Fe/H] = 0.08 ± 0.08 (TTriaud et alJl20Toh . In 
our modeling we included photometric data f r om p re- 
vious publications including iSouthworth et aD (120100 - 
who converted the t iming s fromlCameron et al.l (12007ft . 
iCharbonneau et al.l (|2007l ). Hrudko vTet al.l (j2009T ) into 
the common time st andard BJD(TDB) as outlined by 
Eastman et al.l (12010ft - as well as one transit epoch from 



Sad a et al l (|20i2f ) 



3. APOSTLE OBSERVATIONS OF WASP-2 

The APOSTLE data reduction pipeline is described 
in prior publications, including our observational tech- 
niques, details of photometric extraction, detrending 
of the lig htcurves, and parameter ization of the tran- 
sit model (jKundurthv et al.ll2013bP ). Specifically to the 
WASP-2 system, we acquired 7 Cousins /-band and 3 
r-band transit sequences between July 2007 and Octo- 
ber 2010. The /-band images were fringe-corrected us - 
ing the techniques described in lKundurthy et all (|2013af ). 
Observations were taken with a variety of instrumental 
settings, staring with 0.5s observations on the nights of 
2007-07-24 and 2007-07-26, and moving to longer defo- 
cussed exposures starting in 2010. For the analysis here, 
all data were binned to equivalent -45s exposures. For all 
observations we used the R = 11.4 comparison star TYC 
522-780-1. 

4. MODEL FITTING 




0.05 0.00 0.05 0.10 
BJD - Mid Transit Time 



Figure 1. Seven I- band and three r-band detrended lightcurves 
of the WASP-2 system. The vertical axis is in normalized flux ratio 
units. The horizontal axis shows time from the mid— transit time in 
days, computed by subtracting the appropriate mid— transit time 
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Reduced lightcurves were detrended and mod- 
eled using the Mult iTransitQuick (MTQ) package 
(|Kundurthv et all I2013H) . We used MTQ in two modes: 
to fit for similar transit depths for data taken in a 
given filter (Multi-Filter model); and to fit each transit 
depth individually (Multi-Depth). The set of parame- 
ters used for Multi-Filter version of MTQ is #Muiti-Fiitcr = 
{t T ,t G ,Dj... NF , vij... NF ,v 2 ,j...N F ,Ti... NT }, where t T is 
transit duration and the limb-crossing duration. The 
per-filter fit parameters, up to the number of filters 
Np, include the transit depth D, and limb darkening 
parameters v\ and v%. Finally, the parameter set in- 
cludes mid-transit times Ti up to the number of transits 
Nt- For the Multi-Depth version we used ^Multi-Depth = 
{t T ,t G ,Di... NT ,vij... NF ,V2j...N F ,Ti... NT }, the main 
difference being we fit for each transit depth separately 
(Nt) instead of per-fi lter (Np). 

As outlined in IKundurthv et al.l (|2013bft and 
IKundurthv et al.l (|2013aft . we did not fit for the limb 
darkening coefficients as they result in poorly converged 
Markov chains. Inst ead we kept them f i xed at values de- 
termined using the iClaret fc Bloemenl (|2011[ ) quadratic 
limb darkening models: u\j = 0.3926 , uu = 0.2166, 
Mi T r = 0.5541, M 2 ,r = 0.1594 (see iCsizmadia et al.1 
120131 for cautions regarding this procedure). Our limb 
darkening terms v i , v 2 a re linear combinations of the 
IClaret fc Bloemenl (|201lD quadratic terms, v\ = u\ + u 2 
and V2 — ui — 1*2. 

Figure Q] presents our detrended data, offset for clar- 
ity, with the best-fit model lightcurves overplotted. The 
root-mean-square scatter about the model fits ranged 
from 472 ppm (2008-10-03) to 2252 ppm (2007-07-26), 
with a median of 626 ppm for the /-band data, and 1146 
for the r-band data. 

High resolu ti on i maging of the system by 
IDaemgen et al.l (|2009f) reveals a faint companion 
within the wings of the WASP-2 host star, which will 
affect our conversion from D to R p /R*. This required 
that we estimate the apparent /-band and r-band 
magnitudes of WASP-2A and the contaminating star 
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(named here WASP-2/C). 

We converted the re ported i 1 and z' magnitudes from 
iDaemgen et al.l (120091) using t he transformatio ns pre - 
sented in iRodgers et all (|2006t ) and I Jordi etall (|2006[ ). 
as well as on the SDSS DR7 webpageB The final esti- 
mates yielded r 2A = 11.68 ± 0.11, I 2A = 11.00 ± 0.11, 
r 2/c = 17.05 ± 0.20, 1 2 /c = 15.39 ± 0.20. The two domi- 
nant uncertainties in thes e conversions are : the a pparent 
brightnesses reported by IDaemgen et al.1 <|2009f ) , which 
are uncertain to 0.1 magnitudes; and for WASP-2/C its 
redder color, which leads to a larger uncertainty in each 
color-term. This result indicated that 0.7% of the stel- 
lar flux in the r-band comes from WASP-2/C, and 1.7% 
in the /-band; the WASP-2b transit depth D increased 
proportionally. The change in the J-band depth was 
approximately four times the parameter uncertainty de- 
termined below, making this a necessary correction. The 
derived values of R p /R± increased by 0.3% and 0.8% in 
the r-band and /-band, respectively. 

For each parameter set (^Multi-Filter and ^Multi-Depth): 
we ran two MCMC chains, each having 2 x 10 6 steps. 
These were cropped at the beginning of the chains, 
where the step acce ptance rate is lower than optimal 
(|Gelman et al.l 120031) . yielding approximately 1.8 x 10 6 
steps per chain used in the subsequent analysis. These 
chains were compared against each other to evaluate the 
Gelman-Rubin R-static (|Gelman fe R ubin 1993) and as- 
sure that the chains sufficiently sampled model space. 

Evaluation of our Multi-Depth fit indicated that the 
depths in the /-band were consistent at Dj = 0.0174 ± 
0.0003 magnitudes, while in the r-band the depth were 
also consistent (D r = 0.017 ± 0.001) but with a larger 
RMS due to having only 2 completely sampled transit^. 
We found no evidence for transit depth variations within 
these data, and present the Multi-Filter fits as our final 
results. 

Fitted-for and derived Multi-Filter parameters are 
presented in Table [TJ with the joint -probability distri- 
butions presented in Figure [5] for the fitted and derived 
parameters. Derived parameters include: R p /R+, the ra- 
dius of the planet in units of the host stellar radius; a//?*, 
the normalized semi-major axis a of the planetary orbit; 
the stellar density p*; planet impact parameter 6; orbital 
inclination i; and orbital period P. We found values of 
the Gelman-Rubin R-static within 10 -3 of 1.0 for all fit- 
ted parameters, indicating sufficient coverage of the chain 
over parameter space. The shortest effective chain length 
is for the time of transit on 2007-07-26 with a length 
of 9419. This night had some of the largest photomet- 
ric uncertainties, and largest overall uncertainty on the 
time of transit. All other parameters have effective chain 
lengths larger than 10 4 , i ndicative of sufficient mixing in 
the MCMC sample (e.g. iTegmark et al.ll200l . Finally, 
the x 2 of the best model fit to the data is 2822.93 for 2809 
degrees of freedom (reduced \ 2 of 1.005), indicating that 
our data (and therefore parameter) uncertainties are well 
understood. 

5. SYSTEM PARAMETERS 

3 http: / /www. sdss.org/dr7/algorithms/jcg_photometric_eq_drl. html 

4 Note that these depths were determined before the corrections 
for WASP-2/C were applied, and thus differ from the final Multi- 
Filter results presented in Table [T] 



Fitted-for and derived system parameters are pre- 
sented in Table [1] along with their uncertainties. We 
note that the reported transit depths are for the WASP- 
2 system only, after correction for WA SP-2/C. For com- 
pariso n, we present the results from ISouthworth et ail 
<|201Ct S10), who observed in Cousins .R-band. Figured] 
presents the joint probability distributions of fitted-for 
(left) and derived parameters (right). The left panel 
shows a much smaller level of correlation between pa- 
rameters, making it an effective basis in which to per- 
form MCMC. Final MTQ times of transits for each night 
of observation are presented in the left-hand columns of 
Tabled 

To model the effects of correlated lightcurve noise on 
our t ransit times, we us ed the Transit Analysis Package 
(TAP; IGazak et alll2012ri. which is an impleme ntation of 
the red-noise model of lCarter fe Wmr] (|2009h . Impor- 
tantly, TAP models the amplitude of "white" (random) 
and "red" (correlated) noise on each night of observa- 
tions. We applied TAP using the period derived from 
MTQ, fixed limb darkening values, and zero eccentricity 
and argument of periastron. For APOSTLE observations 
of WASP-2, the median red (white) noise contributions 
were 0.002 (0.0006) magnitudes in the /-band, and 0.006 
(0.0006) magnitudes in the r-band. This indicated that 
while the statistical noise in the two datasets is compa- 
rable, the correlated noise dominates, and is larger in 
the r-band data. This suggests that the correlated noise 
comes from intrinsic stellar variability (either of the tar- 
get or comparison star), or perhaps that the nights of j — 
band observation happened to have larger t ime- variation 
in atm ospheric molecular water absorption ( Stubbs et al.l 
2007). The TAP system parameters are included in Ta- 
ble [I] and times of transit in the right-hand columns of 
Table [U While the TAP uncertainties are typically larger 
than those from MTQ, in all cases the APOSTLE analysis 
yields an improvement in precision over previous mea- 
surements. 

5.1. Transit Timing Analysis 

Using the above analysis, we found a revised ephemeris 
for WASP-2b of 

P = 2.152220976 ± 0.000000305 days 
TO = 2453991. 5148944± 0.0001232BJD 

using APOSTLE results combined with Sad a et all 
(12012ft and the non- amateur results presented in 
Southwor th et al.l (|2010ft . for 17 epochs overall. As Fig- 
ure [3] indicates, there is large scatter in the observed mi- 
nus computed transit times (O-C diagram), with a re- 
duced x 2 of 4.7 for the TAP results (7.0 for MTQ), but it 
is difficult to claim a detection of coherent transit timing 
variations due to the sparse sampling. 

To quantify the significance of this signal, we per- 
formed a 3-parameter sinusoidal fit (period in days, am- 
plitude in seconds, and phase offset as a nuisance param- 
eter) to the O-C data for both the TAP (and MTQ) timings. 
This yielded a A% 2 improvement of 28.2 (47.7), with 
amplitudes of 32 (34) seconds and periods of 389 (437) 
days. We next performed 10 5 random reassignments of 
the timing data to the epochs of observations (i.e. we 
kept the APOSTLE time sampling of any putative TTV 
signal, but shuffled the observed amplitudes). The fitter 
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Figure 2. Plots of the joint probability distributions (JPD) of parameters from the Multi-Filter chains with fixed limb darkening. The 
left panel shows the fitted-for parameters, which are weakly correlated. The right panels shows the derived parameters, which tend to show 
larger correlations. Table [T] provides the relevant units for each parameter set. 



Table 1 

WASP-2 System Parameters 



Parameter 


Value 


TAP 


S10 


Unit 




MTQ 6 Multi 


-Filter Parameters 






ta 


0.0161±0.0002 






days 


tr 


0.0591±0.0001 






days 


D (D 


0.0178±0.0001 








n mq4+ 0002 

u.ui»4_ 0003 








vi(J) 
v 2 (T) 


(0.6092) 
(0.1760) 








ui(r') 
«a(r') 


(0.7135) 
(0.3947) 









Derived Parameters 



(Rp/R*)(i) 
(Rv/R*)(r) 

a/R* 

P* 
b 
i 

P (2.1522 days +) 



0.1315T0.0003 
0.1362^ 



h0.0007 
-0.0009 



8.06T0.04 
2.14T0.03 
0.719T0.003 
84.89±0.05 
1812±26 



0.1317T0.0004 
0.1359±0.001 

7.99±0.06 
2.08±0.05 
0.723±0.005 
84.81±0.08 



0.1326 ±0.0007 
8.08 ±0.12 
2.15±0.09 

84.81 ±0.17 
1852±34 



g/cc 

°(deg) 
milli-sec 



was initialized to the period corresponding to the peak 
of the shu ffled-TTV periodogram, determ ined using the 
method of iZechmeister fc Kiirsterl ([2009). We then ex- 
amined what fraction of these random assignments al- 
lowed a A% 2 equal to or larger than that observed. This 
provided an estimate of the false alarm probability for 
any potential transit timing modulation. We found that 
42% (9.7%) of the random shuffles yielded A\ 2 improve- 
ments at an amplitude equal to or larger than that ob- 
served. While the MTQ results were significant at the 1.7 a 



level, the TAP results are more realistic given the cor- 
related noise in our data. Thus, these data provided 
only marginal evidence for transit timing variations in 
the WASP-2 system. 

6. RESULTS AND CONCLUSIONS 

We presented observations and analysis of 10 transit 
timings of the WASP-2 system by the APOSTLE pro- 
gram. After an extensive treatment of the data to ad- 
dress and model the intrinsic systematic errors, we un- 
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Table 2 



APOSTLE Transit Times for WASP2 


Epoch 


TO (MTQ) 




TO (TAP) 


TO 




2,400,000+ (BJD) 


(BJD) 


2,400,000+ (BJD) 


(BJD) 


146 


54305.73863 


0.00035 


54305.73862 


0.00019 


147 


54307.89212 


0.00063 


54307.89200 


0.00065 


185 


54389.67652 


0.00018 


54389.67646 


0.00024 


298 


54632.87724 


0.00010 


54632.87720 


0.00019 


336 


54714.66135 


0.00014 


54714.66158 


0.00018 


349 


54742.64007 


0.00006 


54742.64011 


0.00009 


362 


54770.61909 


0.00008 


54770.61897 


0.00016 


646 


55381.84717 


0.00092 


55381.84757 


0.00072 


678 


55450.72118 


0.00024 


55450.72108 


0.00029 


691 


55478.70099 


0.00025 


55478.70063 


0.00044 



there is not incontrovertib le evidence of coherent tran- 
sit timing librations (e .g. lEibe et all l2012t iSada et all 
l2(Hl IHover "eTaLl l2L)ll . A primary reason for this is 



Transit Timing Residuals for WASP-2 



Charbonneau (2007) 
APOSTLE (TMCMC) 
Hrudko™. (2009) 
Southwcrth et al (2010) 
Sadaet ale. (2012) 
Collier Cameron (2007) 
APOSTLE (TAP) 



fit 



*5f 



100 200 300 400 500 600 

Epoch 

Figure 3. The observed minus computed transit times for 
WASP-2b. Values from APOSTLEs MTQ analysis, TAP, and pre- 
vious literature are plotted. The horizontal axis represents the 
transit Epoch. The zero-line ephemeris is described in Section |5. II 



dertook a MCMC modeling analysis to understand the 
parameter uncertainties and correlations. As in previ- 
ous publications, a model analysis incorporating tran- 
sit duration tx, limb-crossing duration tg, and tran- 
sit depth D yielded weakly correlated parameters (Fig- 
ure [5]). We corrected for a previously reported object 
withi n the photometric ap erture of the WASP-2b host 
star (jDaemgen et al.l [200 9) to yield the system param- 
eters reported in Table [TJ The uncertainties on these 
paramet ers are smaller than tho se reported in previous 
studies (jSouthworth et all [2010). We disentangled the 
random an d correlated noise in each lightcurve using the 
method of iGazak et al.l (|2012D . and found that the red 
noise dominates scatter in the lightcurves, at a larger 
amplitude in the r-band than in the J-band. 

The depths of transit coming from this analysis did not 
show significant time dependence. However, the times of 
transit shown in Figure [3] show scatter larger than the 
uncertainties. While a basic sinusoidal fit to this signal 
provided a significant improvement in % 2 , we achieved a 
similar goodness of fit in 18% of random reassignments 
of the transit timing variations to the APOSTLE epochs. 
We therefore cannot conclusively report evidence of tran- 
sit timing variations in the WASP-2 system. 

These results mirror many of those reported in the field 
of ground-based transit timing measurements, where the 
O-C diagrams show reduced \ 2 larger than 1.0, but 



sparse sampling of transit epochs, due to weather, day- 
time, and seasonal considerations. In contrast, the con- 
tinuous observations afforded by the Kepler spacecraft 
provide complete sampling of the TTV signal, with both 
high signal-to-noise per transit and a large number of 
transits observed. Resolved Kepler TTVs have sufficient 
overall signal-to-noise to verify and weigh planets in the 
system. Even with tight control of experimental system- 
atics, ground-based observations are intrinsically limited 
by the spatial and temporal stability of the Earth's at- 
mosphere; metrology at the level required to achieve cal- 
ibration at the level of Kepler is currently not available. 

It is likely that the lack of TTVs being detected 
from the ground is due to a fundamental property of 
the systems that are being studied. The vast major- 
ity of ground-based follow-up is focused on hot-Jupiter 
systems, since resolving their transit depths is achiev- 
able even with modest aperture telescopes and non- 
ph otometric obs e rving c onditions. However, as outlined 
by ISteffen et al.l (|2012bT ) , hot- Jupiter systems observed 
by Kepler also fail to show detectable TTVs, or any ev- 
idence for being in a multi-planet system. Their neigh- 
bors in exoplanet parameter space, loosely termed warm- 
Jupiters and hot-Neptunes, do show evidence for both 
TTVs and transiting companions. What this suggests is 
a unique dynamical pathway for the formation of con- 
temp orary hot-Jupiters, such a s multi-planet scattering 
(e.g. [Beaugc fc Nesvornvl I2012D leading to the ejection 
of lesser bodies from the system. It is thus possible 
that those systems that are most easily observable from 
the ground are also those that will fail to exhibit TTVs. 
Given the longevity of ground-based observing resources 
compared to space-based ones, it may benefit the field to 
undertake longer-term observations of systems to resolve 
secula r trends such as those arising from st ellar binar y 
hosts (|Ford et al.ll2000D or tidal orbital decay (|Hutlll981l) . 
The resolution of TTVs from the ground may also require 
observations from sites designed for the minimization of 
aerosols and water vapor, such as Llano de Chajnantor 
in the Atacama desert of Chile, to reduce atmospheric 
red noise and make resolving transits of other classes of 
exoplanet systems feasible. 

While efforts to resolve TTVs from the ground have 
been largely unsuccessful, the drive to follow-up detected 
exoplanet systems has enabled significant leaps in the 
understanding and application of high-precision relative 
photometry. These advances are now extending to the 
broad application of time-domain spectroscopy during 
transits to map out ex oplanet features in absorption (e.g. 
iCrossfield et al.ll2011h . and have overall been a boon to 
the field of time-domain astronomy. 
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